Full counting statistics in a disordered free fermion system 
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The Full Counting Statistics (FCS) is studied for a one-dimensional system of non-interacting 
fermions with and without disorder. For two unbiased L site lattices connected at time t — 0, the 
charge variance increases as the natural logarithm of t, following the universal expression {SN 2 } ~ 
Jj-logi. Since the static charge variance for a length I region is given by (SN 2 ) w ^j-logZ, this 
result reflects the underlying relativistic or conformal invariance and dynamical exponent z = 1 
of the disorder-free lattice. With disorder and strongly localized fermions, we have compared our 
results to a model with a dynamical exponent z / 1, and also a model for entanglement entropy 
based upon dynamical scaling at the Infinite Disorder Fixed Point (IDFP). The latter scaling, which 
predicts (SN 2 ) oc log log t, appears to better describe the charge variance of disordered 1-d fermions. 
When a bias voltage is introduced, the behavior changes dramatically and the charge and variance 
become proportional to (logt) 1 '^ and logi, respectively. The exponent ijj may be related to the 
critical exponent characterizing spatial/energy fluctuations at the IDFP. 
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I. INTRODUCTION 

Full Counting Statistics (FCS) in fermion systems [l[ 
have attracted attention recently for their connection 
with entanglement entropy Specifically, it has been 

shown that the von Neumann entanglement entropy asso- 
ciated with the coupling between two conductors may be 
expressed as a series involving the cumulants of conserved 
charge in one of the conductors @, Q. In some sense, 
noise in an unbiased junction between two conductors 
may be thought of as the propagation of entanglement 
after the junction is closed. In a gapless, disorder-free 
fcrmionic system, it appears that the rate of propagation 
is effectively the Fermi velocity, vp- It is natural, then, 
to ask if entanglement propagates when there is strong 
disorder and the fermions are completely localized? 

Consider a situation where the conductors are 1- 
dimensional and the barrier between them reflectionless. 
In the unbiased case, the second cumulant of charge and 
entanglement entropy are simply proportional (all higher 
cumulants vanish) . For two 1-dimensional noninteracting 
fermion systems of length L, connected at time t = 0, 
the charge variance increases logarithmically in t, follow- 
ing the universal expression (SN 2 ) \ logt, for t much 
shorter than the ballistic time to encounter the boundary, 
tb ~ L [l[. However, the equilibrium charge variance, at 
zero temperature, for a length / region where I << L is 
given by a similar result: (SN 2 ) « -V log/. These results 
may be roughly interpreted in the following way: after 
time, t, charge fluctuations only reflect the participation 
of those fermions within a distance I = vpt of the junction 
Entanglement entropy, 5, of the half-chain exhibits 
the same scaling in space and time but with a different 
universal prefactor: S = ^-(SN 2 ) Q. Whether calculat- 
ing noise or entropy, these result reflects the underlying 
relativistic or conformal invariance and dynamical expo- 
nent z = 1. 

Although the calculations of entropy and noise are 



quite different, both results express underlying univer- 
sality unique to 1-dimensional (critical) systems. It is 
also well known that some features of universality ex- 
tend to entanglement entropy computed in 1-dimensional 
disordered systems. Disorder is strongly relevant in 1- 
dimension and in the case of noninteracting fermions it 
results in complete localization, described by the Infinite 
Disorder Fixed Point (IDFP) ) [5]. Remarkably, entan- 
glement entropy retains its universal logarithmic form 
but with a modified central charge c c fr = c log 2 [6] . Sim- 
ilarly the equilibrium charge variance of a size I region is 
given by (SN 2 ) w ^logZ. 

Since even fully localized fermions exhibit log I entropy 
and noise in equilibrium — similar to delocalized ones — 
one might expect that given enough time the system 
would relax to a state resembling the IDFP with long 
range entanglements straddling the junction after it is 
closed. If nothing else, the fermion wavefunction must 
maintain antisymmetry as fermions cross the junction, 
which is generally a sufficient condition for entanglement. 

In this manuscript we computationally explore the 
temporal development charge fluctuations when strong 
disorder is present. Questions similar to these have 
been addressed in an interesting analytical study of Lieb- 
Robinson type bounds in spatially disordered systems 
Q and in computational works on entanglement after 
a quench || |9|. Very recently, the growth of entan- 
glement after a quench has been studied in connection 
with the role of interactions [lOj. We also note that for 
strictly 1-dimensional conductors, disorder is always rel- 
evant and one would expect the long time limit of FCS to 
be controlled by IDFP physics. For this reason a numer- 
ical analysis of FCS in a disordered 1-dimensional lattice 
might be relevant for experiments on noise in quasi 1- 
dimensional conductors. These results may also be rele- 
vant to recent experiments on cold atom gases in random 
lattices fill. 
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II. THE MODEL 

Consider two one-dimensional noninteracting systems 
of spinless fermions (subsystems A and B), connected by 
a weak link at one site (fig. 1). The model hamiltonian 
is: 

H = £ K xy (c?c% + rfc«) (1) 

(x,y)a=A,B 
I A\ B | B\ A\ 

- w(c 'c +c 'c ) 

= Ho + H w (2) 

where a; and ?/ are one-dimensional site indices and A 
and B denote the two systems with nearest-neighbor hop- 
ping amplitudes n xy and weak link amplitude, w. Each 
subsystem consists of L sites and is taken to have fixed 
boundary conditions at its ends. Furthermore, we will 
restrict our considerations to the case of L fermions in 
2L sites, c x (c" t ,) destroys (creates) fermions at site x in 
subsytem a and obeys the conventional fcrmion algebra. 
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FIG. 1: Two fermion lattices with random hopping ampli- 
tudes K xy , coupled by a weak link with amplitude w. 

Without disorder, all K xy will be set to k = 1. When 
disorder is included in our calculation, the set {k xv } will 
be chosen from a uniform random distribution over the 
interval [k — 5k, k + Sk] . If the {n xy } are chosen indepen- 
dently for the two subsystems, A and B will be at dif- 
ferent chemical potentials. Given exactly L/2 fermions 
occupying each subsystem initially, there will be a net 
particle flux in some direction as the system evolves away 
from its initial condition. To study solely the fluctuation 
in the number of particles in one subsystem, we adopt 
identical realizations of disorder in A and B subsystems. 
Although unphysical, this simplifies the analysis of the 
fluctuations. 

At time t < 0, the two subsystems A and B are discon- 
nected and prepared in the ground state of noninteracting 
fermions: 

L/2 

U') = I[cfc^\Q} (3) 

fe=i 

where (without disorder) the operator ct is related to c x 



in the usual way: 

v fc=l fe=l 

With disorder, (x) will be an eigenfunction of the LxL 
matrix n xy : 

L 

^ K xy (f) k (y) = e k (f) k (x) (5) 
y=i 

It will be useful later to rewrite H in the momentum 
basis 

H = jZe k {c^c A k +c^c B k ) (6) 
fe=i 

k,k' = l 

Ki^ci (?) 

kk' = l;aa'=A : B 

where M kk < = w<j) k (l)4> k ' (1). In the last equality, hf!%, 
may be thought of as unfolded into a 2L-dimension single 
particle hamiltonian matrix. 

At time t = 0, the subsystems are connected with am- 
plitude, w. The probability, foAr(t), that TV particles are 
found in system A at time t is commonly expressed in 
terms of a generating function: 

L 

X (X,t)= 5>aK<K ajv (8) 

N=l 

Following Klich [l2j], the generating function %(A, t) cor- 
responds to the following quantum expectation value, 

X (X,t) = (^e lX ^=i c ^ c *U\iP) (9) 

= {ip\U^e iX ^« a >kk> P M' < >c *>U\ip) (10) 

= {^\e iX T Ja c.'M P uS'^ ca H ca k ''\^) (ii) 

where U = e~ lHt is the full many body evolution oper- 
ator. In the second equality, P^Q is the 2L-dimension 
projection operator, written in momentum space, onto 
the A subsystem. %(A, t) is then expressed (the third 
equality) in terms of the time evolution of the projection 
operator in the single particle Hilbert space (dropping 
indices for clarity): P{t) = e lht p e - lht n 

Now the expectation value with respect to \tp) is 
replaced with a trace over an initial density matrix, 
p = -^e~P H ", corresponding to two independent Fermi 
systems (A and B), at inverse temperature, /3, where 
Z = tre - ^ . The Klich trace formula is then applied in 
the zero temperature limit yielding, 

X (X,t) = lim Itr(e-^V^)) (12) 

p— >oo Z 

= det[l - N + Ne iXP W] (13) 
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where N is the projector onto the Fermi sea of the two 
subsystems; in matrix form, N = 5 aa >6f.k>9(k — L/2). 

Using the property P 2 (t) = P(t), x('M) ma y be ex- 
pressed in terms of the L nonzero eigenvalues of NP(t) 
denoted by {p m {t)}. 



(14) 



The probabilities {&Ar(i)} may be computed from eq. 
(TTJ} using a clever recursive scheme described by Schon- 
hammer in [IH, [l4[ . 
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FIG. 2: (Color online) Natural logarithm of the probability 
distribution 6jv(i) for times t = 1,9, 18 (o, o, n, respectively). 
For t = 18, the probability distribution is very close to a 
Gaussian distribution in agreement with the known behavior 
for a reflectionless barrier as reported in [l[ . 



III. RESULTS FOR FULL COUNTING 
STATISTICS WITHOUT DISORDER 

In the L — > oo, continuum limit of the disorder free 
model described above exact results are known for the 
FCS generating function. When the junction between 
A and B is reflection-less and unbiased, the generating 
function is exactly Gaussian with a variance that depends 
logarithmically on time [l|: 



■ log Kt 



(15) 



Since the probability distribution bpjit) is related to 
x(A,t) by Fourier transform, &/\r(t) follows a Gaussian 
distribution with respect to N for all times. Fig. 2 shows 
results for log 6jv(*) computed on an L = 100 site, ordered 
lattice. For the representative times shown, the distribu- 
tion is Gaussian with a variance that increases in time. 
Noting the vertical scale of the graph, the variances all 
involve fluctuations of 0(1) fermions. 

To examine the Gaussian probability distribution given 
by eq. (|T5"j). we calculate the variance, 



{SN 2 (t)) ee J2 N 2 b N (t) - (J2 Nb N (t)f 



(16) 
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FIG. 3: (Color online) Variance as a function of time, t, for 
an L — 600 lattice. Calculation for variance in translation 
invariant lattice (o) exhibits universal — y logt behavior. Vari- 
ance computed for disordered lattices with disorder strength 
5k — 0.3 (n) and 5n = 0.5 (o) shown for comparison. 



For a system without disorder, we first look in the 
regime 1 << t << kL where the ballistic propagation 
of fermions to the boundary is irrelevant. The topmost 
data in fig. 3 shows the time development of the variance 
in an L — 600 site system for t < 100. A least squares fit 
yields 



(5N 2 (t)) w 0.282 + 0.0985 logi 



(17) 



agreeing within less than 3% of the universal value for 
the prefactor of the logarithmic term: 1/7T 2 « 0.1013. 
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FIG. 4: Variance as a function of time for an L — 100 site 
lattice showing features of conformal invariance. Solid line is 
graph of equation (|18[) . 

When the variance is computed in the regime t £ 
[0, kL] (fig. 4), where the counting statistics are affected 
by the boundaries, the variance reaches a maximum at 
the subsystem half-size, L/2. In terms of the relation- 
ship between noise and entanglement entropy, this result 
might be anticipated from the CFT expression for en- 
tropy after a local quench [l5j]. Klich and Levitov Q 
showed that for a reflectionless barrier, entanglement en- 
tropy at time t after the channel is closed is given simply 
by the second cumulant of the FCS for same geometry. 
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Fig. 4 shows the variance computed for an L 
tice; the data is well fit by an expression, 



(5N 2 (t)) = 0.273 



1 , ,L 7rf 

— log - sin — 



100 lat- 



(18) 



This result agrees with the analytical time dependence 
for entanglement entropy after a local quench in 15| and 
earlier numerical results that anticipated the CFT result 
(second reference of [HI)- We also note that this expres- 
sion may be identified with equation (21) of reference 
by replacing the the spatial interval L with the time 
window, T, as defined in the reference. 



IV. EQUILIBRIUM CHARGE FLUCTUATIONS 

To help frame our results about FCS with disordered 
fermions we will first discuss equilibrium charge fluctu- 
ations. It has been noted that the logarithmic behavior 
of the variance, illustrated in eq. (JT7J) , is a reflection of 
the underlying 1-d nature (in time) of the FCS problem 
in any dimension [l[. Due to relativistic invariance, the 
time scale, t, sets a maximum length scale, I = vpt out- 
side of which number fluctuations are uncoupled. Since 
number fluctuations in a length I segment of a 1-d chan- 
nel are proportional to ^logZ, the charge fluctuations 
in a channel closed for a period of time t cannot exceed 
■^2 log t+c' , where c' is a nonuniversal constant. Equation 
(I18p is also a reflection of the same space/time symmetry 
through the stronger constraint of conformal invariance. 
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FIG. 5: Equilibrium charge fluctuations for a disordered half- 
chain of length I with disorder strength 8n — 1.0 averaged 
over 2000 realizations (o) and disorder-free half chain (n) . 

Although it is a well known result, in Figure 5 (lower 
data set) we compute the charge variance of the A sub- 
system in an L-site disorder free system. The equilibrium 
charge fluctuations are computed from the fermion cor- 
relation function, g(xy) — (c(x)c*(y)) evaluated in the 
ground state of L/2 fermions in an L site system: 



(N 2 



x,y£C 



A least squares fit of the data yields the dependence on 
I, the number of sites in subsystem A, 



(5N 2 )=a+—\ogl 



(20) 



close to the universal result ^2 log^ (a ~ 0.10, b w 0.53). 
The factor of 1/2 appears in front of the logarithm be- 
cause there is only a single boundary. 

For one dimensional noninteracting fermions, disorder 
is a relevant perturbation and leads to the IDFP corre- 
sponding to localized fermions. It is known that entangle- 
ment entropy of a conformally invariant system becomes 
modified with disorder, the entropy following the same 
logarithmic behavior but with a modified central charge 
c c ff = c log 2 [6| . It was shown in [la ] that the variance 
of the z-componcnt of spin in the Random Singlet Phase 
(RSP)[l6| is proportional to the entanglement entropy; 
based on this result, the charge variance for fermions fol- 
lows equation (|2U)) with b = i log 2 « 0.347. Our small 
disordered lattice calculation, averaging over 2000 real- 
izations of disorder (upper data set of figure 5), agrees 
reasonably well with the known result yielding a rs 0.20 
and b ~ 0.365. Since charge fluctuations at the fixed 
point should be independent of the microscopic disorder 
strength, we used a strong disorder parameter 5 k = 1.0; 
presumably our close agreement reflects that we are close 
to the fixed point with this choice. 



V. RESULTS FOR FCS WITH DISORDER 

The IDFP corresponds to completely localized 
fermions. Therefore, one would expect that the heuristic 
argument for the \ log t dependence of the variance — 
specifically, that FCS corresponds to the number fluc- 
tuations in a L — vpt size spatial window — should fail. 
Since the gate time t sets an effective ballistic length scale 
L = Vpt for fermions in a translation invariant lattice, for 
weak disorder one might expect a ballistic regime wherein 
the variance behaves as i log t for vpt < £ where £ is the 
elastic mean free path or localization length. Beyond this 
ballistic regime, at VFt > £, the variance would saturate 
to a value proportional to log £. 

For a system with disorder, we first look in the regime 
1 << t << kL where propagation of fermions to the 
boundary is believed to be irrelevant based upon our 
prior results without disorder. The two lower data sets 
in fig. 3 shows the time development of the variance in 
an L = 600 lattice for two different disorder strengths, 
5k = 0.3, 0.5. Clearly there is a large reduction in the 
variance and a deviation from the universal result for the 
disorder- free system (upper data set). 

The results of a disorder FCS calculation for an L = 
600 lattice with 5k — 1.0) extended to t = 10 4 is shown 
in figure 6. We will discuss the data by fitting it to an 
expression of the form: 



(5N 2 (t))=a- 



■logi 



(21) 
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with 6 = 1 corresponding to the variance without disor- 
der, (k = 1 has been suppressed in the argument of the 
logarithm.) 

In the presence of strong disorder, the charge variance 
appears to be logarithmic in time, but with a prefac- 
tor significantly smaller than that of the disorder-free 
result. A least squares fit to the data of figure 6 gives 
a = 0.27, b = 0.138. 




FIG. 6: FCS computation of charge variance for a disordered 
L = 600 lattice with 8k — 1.0. Results were averaged over 
1000 realizations of disorder. 

Figure 7 shows a comparison of the time dependence of 
the variance for three disorder strengths, 8 k = 0.5,0.75 
and 1.0. Although it is difficult to make a definitive in- 
terpretation, for weak disorder there appears to be a bal- 
listic regime for t < 10. At longer times the variance for 
the weakest disorder strength, 5 k = 0.5, seems to show 
logarithmic behavior with an initial slope corresponding 
to b » 0.394 and a w 0.301. At longer times, t « 100, 
there appears to be a crossover to the logarithmic behav- 
ior associated with the strong disorder of figure 6. 
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FIG. 7: Comparison of variance for an L = 600 lattice with 
three disorder strengths 8k — 0.5 (o), 0.75 (n) and 1.0 (o). 

We first attempt to interpret this result by associating 
a dynamical exponent z with the space and time scaling 
for the charge variance of disordered fermions. Taking 



the universal relation for the equilibrium charge variance 
at the IDFP (b = i log 2 in equation (|2TJ)) , corresponding 
to the single boundary case) as the nominal spatial be- 
havior, we write the temporal charge variance as follows: 



(6N 2 (t)) =a- 



log 2 
2tt 2 



log* 



l/z 



(22) 



If systems A and B are connected at time t = 0, t 1 ^ 
represents the average distance to which fermion entan- 
glements have propagated; z = 2 would correspond to 
diffusive behavior, perhaps indicating diffusive propaga- 
tion of the underlying entangled quanta. Using the fit to 
the data of figure 6, we find z « 2.5. However there is a 
reasonable question about our association single bound- 
ary RSP results for entanglement and noise (see the dis- 
cussion in reference Q)- Comparing equations (fT8")l and 
(|20| we see that in the clean case, dynamic fluctuations 
must be compared to the double boundary static case to 
yield dynamical exponent z = 1. This feature is a conse- 
quence of the relevant conformal map used in the quench 
calculation. With disorder and the breakdown of confor- 
mal invariance, it is not clear that comparison with the 
double boundary case is appropriate. 

Recently, Igloi and coworkers [|[ have proposed a dy- 
namical scaling at the IDFP to explain their results on 
entanglement entropy after a quench. Rather than re- 
placing the spatial scale I with t 1 ^ as in equation (|22|). 
they proposed I ~ (log i) 1 /^. This scaling may be un- 
derstood in terms of the Random Singlet Phase: bonds 
spanning a length I have an exponentially small energy 
scale E ~ exp (— P"), where tp — 1/2. Associating the 
quench time, t, with an inverse energy reproduces the 
proposed scaling. Combining I ~ (logi) 1 /^ with equa- 
tion (f2"0"|) . the variance is then expected to behave as: 



(5N 2 (t))=a- 



1 b , , 
— -log logf 



7T if 



(23) 



Looking at figure 6, the variance appears logarithmic 
over four decades; however extending our calculation out 
to double exponentially long times (as in reference Q) 
and performing more disorder averaging over fewer points 
reveals a deviation from logt behavior. 

Figure 8 shows the variance as a function of log log t 
computed for several lattice sizes ranging from L = 48 
to L = 600. The envelope defined by these data sets is 
approximately linear and saturation of the variance ap- 
pears to occur at doubly exponentially long times. The 
saturation values appear to increase linearly in logL. To 
extract the IDFP exponent, ip, b appearing in equation 
(|23p must be independently determined. Direct analy- 
sis of the variance versus log I for a disordered half-chain 
of length I yielded a value of b ~ 0.365 (equation ([2TJ1) 
and figure 5) close to the analytical result b — log 2/2. 
Following Q we have also determined b from the satura- 
tions values of the variance implied by figure 8. Figure 
9 shows the variance for three largest time values plot- 
ted against logL for the system sizes, L, considered in 
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systematically larger, a result also found in reference [8[. 

We note that there appears to be systematic variations, 
independent of lattice size, in the variance for small times 
and for times larger than the saturation time, t sa t. For 
instance, there is a peak/dip structure at t ~ 4.7 — 6.9 
and another at t ~ 42 — 106 that appears in all the lattice 
sizes studied. There is also a systematic decrease in the 
variance at log log t ~ 3.3 — 3.8, beyond the saturation 
time for each lattice. We speculate that the feature at 
t ~ 4.7 — 6.9 has to do with our choosing identical re- 
alizations of disorder for the A and B subsystems. As 
discussed in section II, this was necessary to equalize the 
Fermi levels in both subsystems. The only reasonable ex- 
planation of the systematic feature at t ~ 42 — 106 is the 
presence of some ballistic fermions. Consistent with ear- 
lier computational work on 1-d localization, we believe 
there may be a substantial fraction of delocalized states 
with localization lengths in the range £ ~ 50 — 100 lattice 
sites (rTj . 
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FIG. 8: (Color online) Variance as a function of 
log log t, for a series of lattice sizes: (from bottom 
to top; number of realization shown in parentheses) 
L = 48 (10000), 64 (10000), 96 (5000), 128 (5000), 192 (3000), 
256 (3000), 600 (3000). In all calculations the disorder 
strength was 5k — 1.0. Dashed line is a visual fit to the 
data envelope: 0.26 + 0.062 log log t. 
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FIG. 9: (Color online) Saturation values of variance. Depen- 
dence of the variance at the three longest time values for the 
lattice sizes from figure 8 upon logL. Dashed line is least 
squares fit: 0.218 + 0.0385 log L 
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FIG. 10: (Color online) Fermions transferred to subsystem 
A (L = 600) as a function of time, t, for several weak link 
amplitudes, w = k, 0. 5k, 0.3k (top to bottom.) Subsystems 
were prepared at t = with an initial bias Nb — Na ~ AN 
'21 The average number of fermions (first cumulant) is seen 
to increase linearly in time up to t — tf> — L, the ballistic time 
to encounter the boundary. For a transparent boundary w = 
K = 1, the first cumulant reaches a maximum (6N) = AiV 
at t — tb = L and then linearly returns to zero at t — 2L. 
Because of partial reflection (w 7^ k), {SN} does not reach its 
maximum until several ballistic time periods, tt, have passed. 



figure 8. Fitting the data in figure 9 to equation (|20|) 
yields a value of b rj 0.380, a value close to that deter- 
mined from figure 5. Finally, iji, from equation (|23p is 
determined to be ijj ss 0.621. Using the value of b from 
figure 5, 1/} w 0.596; using the value of b from the exact 
IDFP (b = log 2/2 w 0.347), V » 0.566. These values of 
if> appear to be close to the analytical value tp = 1/2 but 



VI. FINITE BIAS VOLTAGE 

Since the fermions in the present study are completely 
localized, applying a bias voltage across the junction does 
not result in a constant current. Nonetheless it is inter- 
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esting to determine the rate of charge transfer and the 
resulting shot noise in the 1-d strong insulating regime. 
Based upon our results for an unbiased junction, one 
might expect that the time dependent variance of lo- 
calized fermions described by the IDFP would follow a 
double logarithmic law whether or not there is a bias. In 
contrast, we have found that applying a bias significantly 
changes the dynamic resulting in a time dependent vari- 
ance that is proportional to logt rather than log log t. 



Because of partial reflection (w =^ k), (SN) does not reach 
its maximum until several ballistic time periods, have 
passed. 

To make connections with the Landauer formula for 
conductance (l8j and the Lesovik formula for non- 
equilibrium noise power we denote accumulated 
charge by q{SN). From figure 10 the initial current in 
the w = k case is I = qAN/L. Noting that the density 
of states for the half filled lattice is L/2itk, we arrive at 
the lattice I-V relation: 



5 
4 

^ 3 



♦ ♦ 











200 400 600 800 



t 



FIG. 11: (Color online) time dependent variance for the same 
parameters as figure 10. For w = K = 1, there is only the 
residual universal noise analogous to equation (|18|) describing 
the unbiased case (compare with figure 4). For w — 0.5k and 
w — 0.3k the agreement of our computation with equation 
(|25|) is quite good. (SN 2 ) should be linear in t and achieve a 
maximum value of |w*| 2 (l - K| 2 )A7v. 

Our numerical calculations employ open boundary 
conditions, as opposed to scattering boundary conditions, 
so a closed circuit bias voltage cannot be applied in the 
same way as analytical calculations on FCS. Instead a 
bias voltage is effected by preparing the A and B sides 
with different initial numbers of fermions and watching 
the system relax. As a benchmark we first consider the 
disorder-free case. 

For the results shown in Figure 10, two L = 600 subsys- 
tems were prepared at time t = with different numbers 
of fermions: N A = L/2 - AN/2 and N B = L/2 + AN/2. 
The initial fermion number difference AN was chosen to 
be AN = 24 and the time evolution was studied for three 
different values of the weak link amplitude, w. Figure 10 
shows that the average number of fermions transmitted 
to subsystem A (the first cumulant) increases linearly in 
time up to t = tb = L, the ballistic time to encounter 
the boundary. For a transparent boundary w = k = 1, 
the first cumulant reaches a maximum (SN) — AN at 
t = tb = L and then linearly returns to zero at t = 2L. 



\w*\ 2 ^V 
Zttk 



*V 



(24) 



where 



is the effective transmission coefficient at the 



Fermi level for the weak link and g* is the correspond- 
ing quantized conductance. Referring to figure 10, when 



k the effective transmission coefficient, 



,* 12 



1, 



as expected. When w = 0.5k and w = 0.3k the trans- 



mission coefficients are found to be \w 



* 2 



0.64 and 



\w*\ 2 — 0.3, respectively. The detailed analytical rela- 
tionship between lattice parameters and the transport 
parameters in equations and (|25[) may be found in 
reference [HI]. 

Adapting the Lesovik formula for noise power [l9[ to 
our lattice calculation, we expect: 



(SN 2 



! (1 



Zttk 



(25) 



Figure 11 shows the time dependent variance for the same 
parameters as figure 10. For w = k = 1, the shot noise 
predicted by equation (j2"5)) is zero and there is only the 
residual universal noise analogous to equation (fl"8|) de- 
scribing the unbiased case (compare with figure 4). For 
w = 0.5k and w = 0.3k the agreement of our computa- 
tion with equation (f2"5"]) is close. (SN 2 ) should be linear 



in t and achieve a maximum value of \w*\ (1 ■ 



\w 



2 )AN 



(converting voltage V to equivalent AN) . Using the val- 
ues for w* found from figure 10 for the w = 0.5k and 
w = 0.3k cases, the maximum values of (SN 2 ) are 5.53 
and 5.04, respectively, in close agreement with figure 11. 

We have computed the Full Counting Statistics of dis- 
ordered lattices at finite bias using the same scheme as 
described before in section II. Two identical realizations 
of an L-sitc disordered lattice are prepared with dif- 
ferent numbers of fermions: Na = L/2 — AN/ 2 and 
N B = L/2 + AN/2 and connected at time t = 0. As 
discussed in section II, choosing different realizations of 
disorder for A and B subsystems would introduce a ran- 
dom bias voltage, V r , in effect. To systematically study 
the effect of a bias voltage given the random fluctuations 
of V T , we would have had to bin our multiple realization 
results by the different values of V requiring many more 
realizations to overcome the large statistical fluctuations 
inherent in a small bin size. 

Figures 12 and 13 show (SN) and (SN 2 ) as a func- 
tion of log t for a strongly disordered set of L — 600 lat- 
tices with several different AA^. Saturation in both (SN) 
and (SN 2 ) appears to occur at exponentially long times, 
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FIG. 12: First cumulant (SN) as a function of logt for a 
strongly disordered set of L = 600 lattices with (bottom to 
top) A TV = 4, 8, 12, 16, 20, 24, 28. Each data set represents 
an average over 3000 realizations of disorder (5 k — 1). Sat- 
uration appears at similar exponentially long times as in the 
unbiased case (figure 8). Inset: log- log plot of the data for the 
three largest AJV to determine exponent ip in equation (127[1 , 
Solid line is a visual fit to the linear portion of the data with 
slope V" 1 = 1-25. 
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FIG. 13: (Color online) Variance as a function of logt for the 
same parameters as in figure 12. Data suggests a log log t be- 
havior for log t < log L analogous to the unbiased case shown 
in figures 6 and 8. The variance appears to be proportional to 
log t for large times and saturate at exponentially long times 
consistent with the unbiased case (figure 8). (SN) (figure 12) 
and (SN 2 ) appear to saturate at the same value of AN /A sug- 
gesting that the saturation counting distribution is a poisson 
distribution. 



similar in magnitude to the unbiased case. Using IDFP 
scaling obtained from the unbiased case, we expect the 
saturation time, t sa t ~ expL^". Therefore logt sa t = 53 
or logi sa t = 45 using the two numerically determined 
values of the exponent ip = 0.621 and ip — 0.596 from 
the unbiased variance data. These saturation times are 
in reasonable agreement with figures 12 and 13. Both 
(SN) and (5N 2 ) appear to saturate at the value of AN /4 
suggesting that the saturation counting distribution is a 
Poisson distribution: 

log X (A,i^c«)«^V A -l) (26) 

Using the scaling at the IDFP, one might expect that 
at a time t after the junction is closed, fermions at a 
distance of I ~ (logt) 1 '^ from the boundary are able to 
make a fluctuation bringing them to the boundary. Then 
the number of charges accessible for tunneling would be: 

7 (loa'i) 1 /^ 
(SN) ~ AJV- - AN n [ 8 (27) 

V L (logisat) 1 ^ 

The curvature evident in figure 12 supports a power law 
interpretation. The inset of figure 12 is a log- log plot 
of the data for the three largest voltages in figure 12, 
however a linear fit gives an estimate of ip — 0.8, sys- 
tematically larger than the earlier values determined by 
saturation data. 



The variance (5N 2 ), shown in figure 13, appears to 
have two regimes. At large times running from approxi- 
mately logt = 8 to log/j sat , (SN 2 ) is proportional to logt. 
At short times, logt < logi, is nearly constant and in- 
dependent of voltage. The variance in this regime is very 
close to the unbiased case (see figure 6). Figure 14 shows 
the variance for an L = 300 lattice with a comparable 
range of A/V; the time at which the crossover takes place 
appears to be smaller — possibly proportional to log L. 

Using the saturation value of the variance, AN /A, and 
the linearity in logt we can formulate an empirical law 
for (SN 2 ): ' 



{§N 2 } ^l AN J^L =L i^sL (28) 

X ' 4 log tsat 87TK 



In forming the first equality, we have ignored the short 
time (logt < logL) behavior. In the second equality we 
have used the IDFP scaling for t sat and expressed AJV in 
terms of a voltage, V using the density of states for the 
ordered lattice. We note that the true density of states 
for a 1-d disordered lattice has a complicated singularity 
at the Fermi level and this second equality cannot be 
strictly true for V ->■ [2p||. 
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FIG. 14: (Color online) Variance as a function of logt for the 
same parameters as in figure 13 except L = 300 and AN = 
4, 8, 16 have been changed to to effect similar voltages as figure 
13. The crossover to linear behavior now occurs at smaller 
times, suggesting a crossover at logt ~ log-L. 



VII. CONCLUSION 

In this manuscript we have made a computational 
study of charge fluctuations across a junction of two 
strongly disordered 1-dimensional fermion systems. At 
zero bias, we have analyzed the variance as a function of 
time using a dynamical exponent, and using a recently 
proposed dynamical scaling at the IDFP [8(. By com- 
paring to the known results for static charge fluctuations 



at the IDFP, we have computed the dynamical exponent 
and found z rs 2.5. However upon examining doubly ex- 
ponential time scales, we come to the same conclusion 
as Igloi, et al concerning the dynamics of entanglement 
in reference Q: a region of spatial scale / relaxes with 
a time scale t ~ exp correspondent with the expo- 
nentially small energy scales characteristic of the IDFP. 
This, in turn, leads to the log logt behavior of equation 
(l23l) . We have also concluded that the saturation of the 
variance in a clean system appears to be well described 
by the universal CFT expression developed in references 
15]. 

We have also studied the charge transfer and variance 
with an initial population imbalance, effecting an initial 
bias voltage. For clean systems our results closely fol- 
low Landauer formula for conductivity and Lesovik for- 
mula for noise, modified for our open boundary condi- 
tion lattice geometry. With disorder the fermions are 
localized and charge transfer may be thought of as ac- 
tivated by quantum fluctuations after the connection of 
the two subsystems. Similar to the zero bias case, charge 
and variance are seen to saturate at doubly exponential 
times reflecting the IDFP. However, in contrast with the 
log log t behavior for zero bias, charge and variance relax 
towards their saturation values as (logi) 1 /^ and logt re- 
spectively. It is interesting to note, even with fully local- 
ized fermions, our results for FCS with an applied voltage 
involves the Fermi statistics implicit in the organization 
of the disordered fixed point. 
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was supported by the Department of Energy DE-FG02- 
08ER64623 — Hofstra University Center for Condensed 
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